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Abstract. The Infrared Space Observatory observed the 
field of the 7-ray burst GRB 970508 with the CAM and 
PHT instruments on May 21 and 24, 1997 and with PHT 
in three filters in November 1997. A source at 60 /Ltm (flux 
in May of 66 ± lOmJy) was detected near the position 
of the host galaxy of this 7-ray burst. The source was 
detected again in November 1997, at a marginally lower 
flux (43 ± 13mJy). A Galactic cirrus origin and a stellar 
origin for the emission can be ruled out on the basis of 
the infrared colours. The marginal evidence for variabil- 
ity in the 60 /im flux between May and November is not 
sufficient to warrant interpretation of the source as tran- 
sient fireball emission. However, the infrared colours are 
physically reasonable if attributed to conventional dust 
emission from a single blackbody source. The probability 
of detecting a 60 /im by chance in a PHT beam down to 
a detection limit of 50mJy is ~ 5 x 10~ 3 . If the source 
is at the redshift of the host galaxy of the 7-ray burst 
the fluxes and upper limits at wavelengths from 12 /im 
to 170 /xm indicate it is an ultraluminous infrared galaxy 
(Lj r ~ 2 x 10 12 Lq). The star formation rate is estimated 
to be several hundred solar masses per year, depending sig- 
nificantly on model-dependent parameters. If this source 
is associated with the host galaxy of GRB 970508, pro- 
genitor models which associate GRBs with star-forming 
regions are favoured. 



1. Introduction 

The recent discoveries of fading afterglows to a number 
of 7-ray bursts (GRBs) have been precipitated by the 
accurate and prompt localisation capability of the two 
Wide Field Cameras (WFC) aboard the 'BeppoSAX' X- 
ray satellite (Boella et al. 1997; Piro et al. 1998a). The 
typical ~ 3' radius error circles, obtained in a matter of 
hours after the occurrence of the GRBs, are easily cov- 
ered by ground-based optical and radio telescopes, allow- 
ing rapid and deep follow-up observations. The reduction 
of GRB error circles by BeppoSAX also made it feasi- 
ble to study the content of these regions at far-infrared 
wavelengths with the European Space Agency's Infrared 
Space Observatory^ ISO ( [Kessler et al., 1996| ). The error 
region of GRB 970402 was the first to be rapidly surveyed 



at far- infrared wavelengths (Castro-Tirado et al., 1998a) 



ISO observed this GRB error circle 55 hours, and again 
8 days, after the burst event and detected no new tran- 
sient sources down to a 5 a limit of 0.14mJy at 12 /im 
and 350 mJy at 170 /im. Details of the ISO observations 
of GRBs may be found in Delaney et al. (1999) . 

2. Observations of GRB 970508 



GRB 970508 triggered the GRB Monitor (Frontera et al 



1997) of BeppoSAX at 21:41:47 UT on May 8, 1997 
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( posta et al., 1997 ). The GRB localisation was subse- 
quently improved by the BcppoSAX WFC2 detector to 



a 3' radius (99% confidence) error circle (rleise et al. 
1997|)~|Follow-up observations by the BeppoSAX narrow- 
field instruments (NFI) then detected a new, fading X-ray 
source, consistent with the GRB error circle, but with an 
improved positional accuracy of 50" radius (Piro et al. 
1998b). Within this error circle an optical transient (OT) 
was discovered whose flux increased for the first two days 



and subsequently decayed ( Bond, 1997 ; Djorgovski et al 
1997fj~[rhe decay phase was characterised by a power law 
depen dence w ith time, f(t) oc t~ 1A days ( Castro-Tirado 



et al., |l998b|) . Near-infrared (2.2/im) emission was also 
observed from the OT source between May 13.25 (K s = 
18.2± .2 ma g) and May 20.21 (K s = 19±0.3 mag) flChary 



et al, |1998 ). Spectra of the OT taken with the Keck II 



10 m telescope revealed the presence of [Oil] emission and 
Fell and Mgll absorption lines consistent with a source at 
a redshift of 0.835, providing the first direct measurement 
of the distance to a GRB (Metzger et al. 1997a, 1997b). 
Radio afterglow was also detected from a source at the 
position of the OT (Frail et al., 1997[ [Bremer et al., 1998| ; 



Galama et al., 1998). Observations by HST have shown 
the OT position to be remarkably well-centred with re- 
spect to the host galaxy ( Fruchtcr et al., 1999| ). The host 



is late-type, with the characteristics of a blue compact 
galaxy ( |Sokolov et al., 199£ ). 

We report here on ISO observations of the field of 
GRB 970508 in the weeks following the GRB and on a 
follow-up observation made in November 1997. 

3. ISO Observations 

The GRB 970508 error circle was not initially visible with 
ISO due to orbital constraints but became visible for 
seven short windows lasting from 3000 to 8500 seconds 
between May 20 and May 26, during which time tar- 
get of opportunity observations were made with ISOCAM 
(|Cesarsky et al., 1996[ ) and ISOPHOT ( |Lemke et a"T 



1996). A second window was available from October 25 



to November 21 1997 and follow-up observations were 
then made with ISOPHOT. Details and results from the 



CAM observations are presented elsewhere (Hanlon et al 



199J). The May observations were centred at RA(2000): 
6h53m46.7s, Dec(2000): +79° 16' 02.0", the cen tre of 
the X-ray afterglow error circle (Piro et al., 1998a). The 



November observations were centred on the position of th e 
optical transient associated with the GRB ( Bond, 1997 ). 
Raster maps of the NFI 50" radius error circle were 



obtained w ith the PHT C100 and C200 detectors QKlaas 
et al., 1994|). The PHT C100 camera is a 3x3 pixel Gc: Ga 



array, with a pixel centre-to-centre distance of 46" in both 
directions. The C_60 filter (reference wavelength at 60 /xm, 
width 23.9 /im) and the C_90 filter (reference wavelength 
at 90 /im, width 51.4 /xm) are associated with the C100 
detector. Observations in these filters employed a 5x5 



position raster grid with a 23" step size in both direc- 
tions. The grid was aligned with the orientation of the 
detector array. The size of the area covered was therefore 
230" x 230". This mode of observing was selected to obtain 
a fully sampled map with a high level of redundancy in 
the map centre and to minimise the uncertainties due to 
cirrus confusion. The redundancy was highest in the inner 
46" x 46" where it was 25. At 60 and 90 /im the FWHM of 
the beam profiles (psf convolved with pixel) are 41" and 
47" respectively. 

The PHT C200 camera is a 2x2 pixel stressed Ge: Ga 
array, with a pixel centre-to-centre separation of 92". A 
5x5 grid raster with a 46" step size in both directions 
was used for the C200 observations in the C_160 filter 
band (reference wavelength at 170 /im, width 89.4 /im). 
The area covered was 368" x 368" with a redundancy for 
the inner 91" x 91" of 16. 

4. ISOPHOT Data Processing 

The ISOPHOT data were processed using the software 
package PIA version 7.3 (Gabriel et al., 1997). All stan- 



dard processing steps were included to obtain the flux 
per raster point per detector pixel before map coaddi- 
tion (Laureijs et al., 1998). This flux was taken to be 



the median value of all signals measured during a given 
raster point integration thereby minimising any positive 
flux bias due to the relatively high number of glitches, in- 
duced by ionising particles, in the data stream. The data 
were flatfielded by imposing the same integrated flux for 
all detector pixel scans. This operation yielded a correc- 
tion factor for each pixel without affecting the total flux. 
The maps were constructed by coadding the flatfielded 
fluxes per raster point per detector pixel in a regular im- 
age grid of 23"x23" image pixels which are aligned with 
the orientation of the detector pixels. 

Photometry on the C_60 and C_90 maps was carried 
out by integrating the flux in an area of 5 x 5 image pixels, 
minus the corner pixels, in the region of interest and then 
subtracting the mean background flux derived from the re- 
maining pixels in the map. The measured signal was then 
scaled up to the full PHT beam-size. Since the maps from 
May and November have different orientations and map 
centres, the 60 /im, 90 /im and merged (combined 60 /im 
and 90 /im) maps were reprojected to have the position of 
the OT as map centre and with zero position angle. These 
maps were used to derive photometric fluxes. 

5. Results 

A point source was detected at 60 /im on May 24, con- 
sistent with the OT position, given the uncertainty in the 
position of the peak of the beam profile with respect to the 
geometric shape of a given pixel (Fig. |l|). For a relatively 
low signal to noise detection such as this (see Table 0), 
which is significant in only a few pixels (not all 9), an 
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offset of even a few arcsec can cause the mapping rou- 
tine to shift the source by an entire image pixel. However, 
it is possible that the offset of the source relative to the 
OT position is real and the consequences of this are dis- 
cussed below. Further observations were made with PHT 



D6 h 53 m 00 s 




Fig. 1. PHT C100 map using the C_60 filter, taken on 
May 24, 1997. The OT position (marked with a cross) is 
consistent with the location of the C_60 source (marked 
with a circle) to within the HPBW of the C100 detector. 
The size of the HPBW is shown on the side of the figure, 
along with the grey scale values. 

in November 1997 to verify the detection. The source was 
again detected at 60 /im (albeit at a lower significance level 
than in May) but no significant excess at the position of 
the OT was detected at 90 /im (see Fig. | and Table Q). 
Although the 60 /im flux is lower in November than in 
May, we do not regard this as providing strong evidence 
for variability, since the two fluxes agree to within 2a. The 
offset of the source relative to the OT is again observed in 
the November map. In an attempt to improve the signal 
to noise ratio of the C100 observations from November 
1997, the data from both filters were merged into a single 
map (Fig. |^). The 60 /im and 90 /xm data were separately 
flatficldcd and the corrected fluxes per pixel per detector 
and per filter were coadded on a common grid. The fluxes 
per filter were weighted by their respective uncertainties. 
The merged map has a noise level intermediate between 
the noise in the individual maps and the significance of 
this is discussed below. In addition, the source is detected 
at a signal to noise ratio consistent with that of the 60 /im 
map alone. This suggests that the source is not significant 
at 90 /im. Upper limits (3a), derived from the 170 /im ob- 
servations in May and November, are given in Tabic |l|. 
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Fig. 2. PHT C100 map using the C_60 filter (top) and the 
C.90 filter (bottom), taken on November 11, 1997. The OT 
position is marked with a cross while the position of the 
C_60 source detected in May 1997 is marked with a circle. 
A source is present in the C_60 map which is consistent 
with both the position of the OT and the C_60 source de- 
tected in May 1997. Note: The orientation of these maps 
is different to that in Fig. 1, due to differences in orien- 
tation and map centres between the May and November 
observations. 



UT Start - End 


Instrument 


Filter 


Flux Density 
(mJy) 


May 1997 








24.095 - 24.106 


PHT 


C_60 


66 ± 10 


24.107 - 24.114 


PHT 


C.160 


< 70 


November 1997 








11.021 - 11.040 


PHT 


C_60 


43 ± 13 


11.041 - 11.065 


PHT 


C_90 


19 ±7 


11.066 - 11.070 


PHT 


C.160 


< 130 


Merged C_60 plus C_90 


32 ±9 



Table 1. Log of ISOPHOT observations of GRB 970508. 
The PHT fluxes correspond to the reprojected maps (see 
text for details). 
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Fig. 3. Merged 60 fim and 90 fim map from November 11, 
1997. The OT position is marked with a cross. The local 
maximum near the position of the OT is marked with a 
circle. There are other sources near the edge of the map 
which are brighter but have lower signal to noise due to 
the lower coverage in those areas. 



5.1. Assessment of uncertainties in the maps 

Per raster point, the integration times were 20 s at 60 /im 
in the May observation and 60 s at 60 /im and 90 /im 
in the November observations. Although the integration 
time per raster point was three times longer at 60 /im in 
the second epoch, the noise levels per raster point in the 
two 60 /im maps are comparable (35mJy/beam in May, 
31mJy/beam in November). The relative increase in the 
noise level of the C_60 measurement in November is most 
likely due to the fact that the observation was made at 
the end of a revolution and was therefore subject to an 
increased gli tch rate as the sate llite re-entered the radi- 
ation belts ( Kcssler et al., 1996] ). The noise level in the 
C_90 measurement was 17mJy/beam. Including the re- 
dundancy, the predicted flux level in the map centre is 
about 12 mjy/beam at 60 /im and 6 mJy/beam at 90 /im. 

The detection limit is also determined by the cirrus 
confusion noise (N), which can be estimated from 



iV(mJy) = 1.08 



( A 



VlOO/im 



2.5 



( 



VIMJy/sr 



1.5 



where A is the wavelength, I v is the background surface 
brightness level and a telescope of diameter 60 cm is as- 



sumed (Helou & Beichman, 1990). The validity of the for- 
mula is discussed in Herbstmeier et al. (1998) , and is prob- 
ably better than an order of magnitude. The detection 
limit also depends on the absolute calibration of the back- 
ground surface brightness. The maps indicate background 
surface brightness levels in the range 3-6MJy/sr which 
correspond to cirrus confusion noise fluxes of N(60 /im) ~ 



1.5-4.4 mJy and N(90 /im) ~ 5.6-15.8 mJy. These numbers 
indicate that most of the noise measured at 90 /im (7mJy 
according to Table |l|) is probably due to cirrus structures. 
However, this is not the case at 60 fan. Therefore the noise 
in the merged 60+90 /im map is higher than the noise in 
the 90 fim map, contrary to the predicted weighted un- 
certainty if the noise is purely statistical. The presence of 
cirrus structures at 90 /im therefore imposes a limit on the 
signal to noise which can be achieved. 

6. Discussion 

6.1. Origin of the 60 fim emission 

Blink comparison of the November 60 fim and 90 /im maps 
(Fig. ||) indicate a striking similarity in background struc- 
ture between them. A correlation analysis confirms that 
there is a structure in the maps (correlation coefficient 
0.62) but that the random noise dominates over the intrin- 
sic variations in the background. The surface brightness 
ratio (in MJy/sr) determined from the maps is consistent 
with that expected from cirrus in the Galaxy (190/160 = 
2.03 ± 0.9). It is therefore plausible that the background 
structures in the maps are due to cirrus. If the 90 /im emis- 
sion at the position of the GRB is entirely attributed to 
cirrus, then the maximum cirrus contribution to the flux 
at 60 /im is ~10 mJy. Subtracting this component from the 
observed fluxes in May and November, an averaged source 
flux of 45±18mJy is obtained. On the other hand, if the 
60 /im flux is attributed to cirrus then the flux at 90 /im 
should be ~ 100 mJy which is clearly not the case. We 
can thus rule out a Galactic cirrus origin for the source at 
60 /xm. 

A stellar origin for the 60 /im source can also be ex- 
cluded on the basis that the ISO colours are incompatible 
with stellar spectra over the full range of spectral types 
( [Smith et al., 19871 ; |Hcnden fc Stone, 1998j ). 



We can also exclude with reasonable confidence the 
possibility that the 60 /im source is entirely due to tran- 
sient fireball emission from the GRB on the basis that 
(a) the flux at 60 /im on May 24 is about two orders of 
magnitude greater than the extrapolated radio-optical af- 
terglow spectrum on the same date and (b) the source 
is detected again in November at a consistent flux level 
(to within 2a). If the 60 /im emission followed the t -1 " 1 
decay of the afterglow, its flux should have decayed by a 
factor of roughly 15 by November. The relativistic shock 
scenario accounts quite well for the features of the broad- 
band spectral energy distributions of GRB afterglows and 
a new physical mechanism, beyond the scope of this paper, 
would have to be developed if the extraordinary excess at 
60 fim was to be connected to the fireball. We now argue 
instead that the infrared colours we observe in the source 
are physically reasonable if attributed to conventional dust 
emission from a single blackbody source. 

Modified blackbody spectra of the form v n x B J/ (T), 
were fit to the range of 60/im/12/im colour ratios (R) 
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obtained from observations in May 1997 (Hanlon et al. 
19990 

F60/mi=66mJy/F12/im=0.065mJy =>• R ~ 1000 
F60/im=76mJy/F12/iin=0.035mJy R ~ 2200 
F60/im=56mJy/F12/im=0.095mJy ^ i? - 600 

No colour corrections were incorporated in view of the er- 
rors on the flux determinations. Estimated temperatures 
for modified blackbody spectra with < n < 2 were then 
obtained for these R values and predicted F90 /im/F60 fim 
ratios were calculated for each derived dust temperature 
and assumed emissivity (Table 2). The temperatures for 



R 


T dust (K) 


Predicted F90 /im/F60 jum 




n = 


n = 1 


n = 2 


n = 


n = 1 


n = 2 


600 


91 


81 


70 


0.79 


0.58 


0.44 


1000 


87 


76 


68 


0.83 


0.61 


0.45 


2200 


82 


72 


65 


0.86 


0.64 


0.47 



Table 2. Modified blackbody spectral fits for different val- 
ues of R, yielding a range of estimated dust temperatures 
and predicted F90 /im/F60 /im colours. The R values were 
derived from the May 1997 CAM and PHT data. 



each dust model are well constrained. Moreover, the pre- 
dicted 90 /im fluxes do not vary significantly for a given 
n even when R changes by almost a factor of 4. These 
results also show that subtracting a lOmJy cirrus com- 
ponent from the 60 /im flux does not significantly affect 
the conclusion. Since no observations were made at 90 /im 
in May 1997, due to orbital constraints, we cannot com- 
pare these predicted ratios with observations at the same 
epoch. However, the fit results can be used to predict the 
90 /im flux in November 1997 by appropriate scaling of the 
60 /im flux ratios from May (66±10mJy) and November 
(43±13mJy). Taking into account the photometric errors, 
the ratio of the May to November 60 /im flux is in the 
range 76/30 (2.5) to 56/56 (1). Multiplying the predicted 
F90 /im/F60 /im ratios in Table 2 by these scaling factors 
we see that the predicted 90 /im flux for November 1997 is 
consistent with the observed 90 /im value for all dust mod- 
els when the scaling factor is 2.5 and for the n = 2 case 
when the scaling factor is 1. Therefore it is not necessary 
to speculate that the infrared source has varied between 
the two epochs because the ISO observations are consis- 
tent with dust emission from a single blackbody with a 
temperature of ~ 70 K and n = 2. The predicted 170 /im 
fluxes ^re always within the upper limits obtained (Ta- 
ble [l]), being a factor of two lower than the 90 /im values. 



and qo=0.5. If the upper limits across the spectral range 
are included and the area under the SED is integrated us- 
ing a power law interpolation between spectral points, an 
upper limit on the bolometric luminosity (0.4-1000 /im) 
of 2.9 x 10 12 L Q is obtained. A value of 2.0 x 1O 12 L 
is obtained if we consider only the ISO data points and 
upper limits as contributing towards the luminosity and 
clearly the emission is dominated by the infrared com- 
ponent. Alternatively, the 60 /im source may be associ- 



Restframe wavelength (/^.m) 



10" 10 



-,28 





1 7Q/j,m 

•4\t 




1350/iftT 
86. B GHz -J/ 


i\ 

90^m '. 




1.49 GHz' 












Arp220 




HK s 



,10 , 11 10 12 , 13 10 1 
Restfrarne frequency (Hz) 



,•5 



Fig. 4. Spectral energy distribution of the 60/im source, 
assuming it is the host galaxy of GRB 970508. Data and 
upper limits are from Bloom et al. 1998; Pian et al. 1998; 
Smith et al. 1999; Galama et al. 1998; Shepherd et al. 1998 
and the present work. The 12 /im upper limit was obtained 
by co-addition of the 2 CAM LW10 May observations. The 
C_60 data point represents the averaged flux from May 
and November. The SED of the ULIG Arp 220 is shown 
for comparison purposes only. 

ated with the foreground absorbing system at z=0.767 in 
which case the infrared luminosity is ~ 1.7 x 10 12 Lq. In 
either case the source is an ultraluminous infrared galaxy 
(ULIG, L ir > 10 12 L Q ) ( |Sanders fc Mirabel, 1996| ). Al- 



though ULIGs were discovered by the IRAS satellite in 
our local universe (z < 0.2) (Soifer et al., 1986) a recent 
re-analysis of the IRAS data has shown the existence of 
ULIGs out to z=0.529 ( |van der Werf et al., 199S| ). The 
extremely rare hyperluminous infrared galaxies (HyLIG 
Li r > 10 13 Lm) have been detected to z=2.2 ( |Rowan- 



If the 60 /im source is associated with the GRB host 
galaxy, at a redshift of 0.835, limits can be placed on 
its bolometric luminosity. A broadband spectral energy 
distribution (SED) for the 60 /zm source, assuming it is 
the ho st galaxy of the GRB, is shown in Fig. fl|. The as- 
sumed cosmological parameters are Hp— 50 km s -1 Mpc -1 



Robinson et al., 199l| ; [Cutri et al., 1994| ). 

There arc two galaxies within 5" of GRB 970508, either 
of which may be responsible for the absorption system 
seen at z=0.767 ( [Pian et al., 1998j ). Due to the large PHT 
pixel size, it is possible that one of these galaxies, rather 
than the GRB host, is the counterpart of the 60 urn source 



The g alaxy to the north-west of the GRB, 'G2' ( Pian ct al 



1998 ) has the colours of a late- type spiral at a redshift of 
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0.7, making it unlikely to be the counterpart, since IRAS 
detected no late-type spirals with L; r > 10 11 Lq ( [Banders 
& Mirabel, 1996). Galaxy 'Gl' to the north-east, has ex- 



is observed at visible frequencies represents emission from 
stars forming near the edges of clouds which can escape 



directly without being reprocessed by the dust (Rowan- 



tremely blue colours indicative of a rapidly star-forming Robinson ct al., 1997). The SFR can be deduced from the 



system (Pian et al., 199S) and we cannot rule out that it 
is the counterpart of the 60 /im source. 

The prototypical ULIG is Arp 220 with a far-infrared 
luminosity > 10 12 L Q . The SED of Arp 220 (Fig. @) con- 
sists of a cold component, corresponding to a modified 
blackbody temperature of ~ 50 K and a warm component 



far-infrared luminosity using the following expressions for 
M*, the rate of star formation per year in units of Mq: 
M, = 2.60/e x £so x 10 -io arK } M* = 9.36/ e x ^ x 10~ 10 

where <j) incorporates a correction factor from a Salpeter 
initial mass function (IMF) to the true IMF and a correc- 



tion if the starburst is only forming massive stars ( Rowan- 



(Klaas til al., 1997). Il has 1 been suggested LhaL Lhc warm 



with ajgggro of 120K poking 5Wwd of 60/llli Robinson ct g~ 1997 fr, e is the fraction ^ 1} £ 



component (L warm = 10 11 Lq) is heated by an active nu- 
cleus, while the cold component (L eo id = 10 12 Lq) is star- 
burst powered ( [Rowan- Robinson fe Efstathiou, 1993| ). 

6. 2. Probability of a chance superposition 

The a posteriori probability of finding a 60 /im source 
within a PHT 41" diameter beam can be estimated by ex- 
trapolating recently determined ISO survey source counts 
at 15 /im, 90 /im and 170 /im. At 15 /im the integral source 
counts at 50mJy from the ELAIS survey are ~ 10 4 per 
steradian while at 90 /im there are 10 5 to 10 6 sources per 



steradi^ii to Lhc same flux limiL ( rtnwaii-Robiiisoii cL al 



199S| ). In a deep survey of the Lockiiian Hole 



source 



counts at 95 /im and 170 /im at 150 mJy extrapolate into 
the same range as ELAIS at 50mJy ( |Kawara et al., 1999| ). 



The consistency of the source counts from yu to 1 ( b /im 
and the increasing source counts from lb/im to VTSjEm 
indicate that an estimate of ~ 10 5 sources per stera- 
dian at 60 /im is reasonable. This implies a probability 
of ~ 5 x 10~ 3 of finding a 60 /im source down to a de- 
tection limit of 50 mJy by chance in a PHT beam of 41" 
diameter. 

Alternatively, modelled source counts at 60 /im can 
be used to derive the probability that the 60 /im source 



is a chance superposition with the GRB host (Pearson 



& Rowan- Robinson, 1996). Taking into account source 
counts due to normal, starburst, Seyfert and HyLIGs, the 
probability of finding a 60 /im source with the observed 
flux within a PHT beam is also ~ 5 x 10 -3 . These esti- 
mates are sufficiently small to hypothesise a physical con- 
nection between the 60 /im source and the host galaxy of 
the GRB. 

6.3. The star formation rate of the 60 (im source 

It has long been recognised that some starbursts are ob- 
scured by dust ( Kcnnicutt, 1998j ) and recent results from 
ISO have shown that such obscuration is widespread and 
important for the history of star formation in the universe 



(Elbaz et al., 1999). Hence far-infrared luminosities pro- 



vide a more reliable estimate of star formation rates in 
galaxies. If the source is a ULIG, then most of its massive 
star formation occurs in dense molecular clouds and what 



cal and ultraviolet energy emitted in a starburst (last- 
ing 1 Gyr) which is absorbed by dust and re-emitted in 
the far infrared. Lgo and L15 are the restframe 60 /im and 
15/xm luminosities {vh v expressed in units of Lq). These 
two measures yield values for M* of between 1900 and 
2200 Mq /year. An alternative estimate which uses the 
infrared (8-1000 /im) luminosity and which assumes star- 
bursts lasting < 1 s years, yields a SFR of ~ 500 MQ/year 
( [Kcnnicutt, 1998| ). 

The identification of a ULIG as the possible parent of 
the GRB host galaxy therefore favours models of GRB 
formation involving compact stellar remnant progenitors, 



such as supernovae (Berczinsky et al., 1996; Woosley. 



1993) or hypernovae (Paczynski, 1998), which place GRBs 



in or near star-forming regions. The possible detection of 
an iron line in the X-ray spectrum of GRB 970508 pro- 



vides independent support for this scenario ( Piro et al 



1999 ). However, the HST data show that the optical tran- 



sient is well-centred on th e host galaxy, which is quite 
blue ( Fruchtcr ct al., 1999 ). Since ULIGs tend to be dis- 
turbed, dusty systems, with a significant fraction under- 
going mergers, the host galaxy would be expected to be 
redder in colour and its morphology more disturbed than 
is observed ( |Clements ct al., 1996 ). Another possibility 
is that the GRB host has been formed by the ULIG, as 
the product of a merger. Small galaxies formed along the 
tidal tails in mergers are likely to detach to become dwarf 
systems with characteristics (e.g. colour, absolute magni- 
tude) similar to the GRB host ( [ganders fc Mirabel, 1996 



[Hunsberger et al., 1996 ). We can speculate in that case 



that the offset observed in the May and November 60 /im 
maps between the 60 /im peak and the OT is genuine, 
arising from the observation that the most luminous tidal 
dwarfs are those at the largest projected distances from 



the parent nucleus (Hunsberger et al., 1996). The galaxy 
Gl may also have been formed in the merger but a direct 
redshift determination would be required to support this 
hypothesis. 

7. Conclusions 

The ISO observations of GRB 970508 place unique lim- 
its on the level of far-infrared emission in the weeks and 
months following the burst event. Non-transient emission 
observed at 60 /im indicates the presence of a ULIG which 



L. Hanlon et al.: ISO detection of a 60 /im source near GRB 970508 



7 



may be the parent of the host galaxy of the GRB. This 
result may have important implications for GRB progen- 
itor models, favouring those which place GRBs in or near 
star-forming regions. 
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